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Abstract. The causal character of the zero-areal-radius (R — 0) singularity 
in spherically symmetric spacetimes is studied. By using the techniques of 
the qualitative behaviour of dynamic systems, we are able to present the most 
comprehensive scheme so far to try to find out their causal characterization 
taking into account, and analyzing, the possible limitations of the approach. 
We show that, with this approach, the knowledge of the scalar invariant m = 
R(l — g^ iU d fl Rd u R)/2 suffices to characterize the singularity. We apply our 
results to the study of the outcome of Black Hole evaporation and show different 
possibilities. In this way, we find that a persistent naked singularity could develop 
in the final stages of the evaporation and we show its distinctive features. Likewise, 
we study the options for the generation of naked singularities in the collapse of an 
object (such as a star) as a means of violating the cosmic censorship conjecture. 



PACS numbers: 04.20.Gz, 04.20.Cv, 04.40.-b, 04.90. +c 
1. Introduction 

Singularities are not part of the spacetime since they are related to diverging 
curvature invariants, to the incompleteness of curves and/or to the lack of tangent 
vectors. At most, it would seem reasonable to say that singularities are situated in 
the boundary of the spacetime, provided a suitable definition of boundary is given. 
In fact, one such definition was first introduced by Penrose in 1963 pp. His idea was 
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to embed the spacetime under study with metric g into another Lorentzian manifold 
(the unphysical spacetime) with metric g conformally, 

g = ^ 2 g, 

so that the causal properties are trivially kept. In this way, the boundary acquires 
causal properties itself which are obtained by its mere examination in the unphysical 
spacetime. Specifically, it becomes now meaningful to give a singularity attributes 
such as spacelike, timelike or lightlike. Furthermore, in spherically symmetric 
spacetimes, where the SO (3) group orbits form a spacelike two surface (the 2- 
spheres), it is possible to perform just the conformal compactification of the 
two-dimensional surface orthogonal to the 2-spheres retaining all the important 
information. This is so because, by means of a coordinate change, the induced 
Lorentzian metric or first fundamental form of the two-dimensional surface can 
always be brought into a conformally flat form 

ds\ d = fi 2 (x , Xi){— dx\ + dx\). 

In this way, it can be naturally embedded in an unphysical two-dimensional 
Minkowskian spacetime (see, for instance, [2]). This allows to draw simple two- 
dimensional diagrams, called Penrose diagrams, from which one can find out the 
properties of the boundary at a glance. In case the boundary is a C 1 curve at a 
given point in the Minkowskian spacetime one could causally characterize the point 
according to the tangent vector to the curve as usual. However, as the cases analyzed 
in the literature show[3], one usually finds boundaries that are just piecewise C 1 and 
which provide us with "piecewise causal characterizations" . 

In order to find out the main features of the local causal character of a 
singular boundary in a spherically symmetric spacetime it is not necessary to follow 
completely the conformal procedure explained above. This is interesting because an 
analytic conformal compactification can only be found for certain particular cases. 
The alternative procedure is based on the general fact that the concept of a null 
geodesic is a conformally invariant one, so that for every null geodesic in the physical 
two dimensional surface there is a corresponding null geodesic in the unphysical two 
dimensional minkowskian spacetime jl]. It seems reasonable, as we will try to show, 
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that the study of the behaviour of the null geodesies in the physical spacetime in the 
neighborhood of a singularity will provide us with information about the behaviour 
of the geodesies around its corresponding conformal boundary and, as a consequence, 
on its causal character. 

On the other hand, we will mainly deal with probably the most interesting 
type of singularities in spherically symmetric spacetimes: The zero-areal-radius^ 
scalar curvature singularities j5]. In order to define the concept we will use the areal 
radius R such that the area of a 2-sphere is 4ttR 2 . Then, we say that there is a 
zero-areal-radius scalar curvature singularity at a point p in R = if any scalar 
invariant polynomial in the Riemann tensor diverges when approaching it along any 
incomplete curve. 

The difficulties to apply the above procedure will lie not in the identification 
of this type of singularities, but in the study of the behaviour of the radial null 
geodesies. In order to carry out this study we will analyze the system of differential 
equations that describe the null geodesies by means of the standard qualitative theory 
of dynamic systems [6j [7J El IH [10] . (Note also that we have written appendices A and 
B with the main results on this subject for the less usual cases). The information 
supplied by this theory about the behaviour of the radial null geodesies will allow 
us to mathematically classify the neighborhood of p in the cases where the theory 
can be applied. As we will see, this together with the correct interpretation of the 
results for the R > region, will provide us with a piecewise characterization of the 
singularity around the chosen point. 

But before putting this plan into practice, let us comment on the previous 

results and the relevance of this analysis. In fact, the study of R = singularities 

has been carried out for important particular solutions. In some cases the conformal 

boundaries have been obtained by using a conformal compactification (for example, 

Schwarzschild's and Reissner-Nordstrom's solutions -see, for instance, [11]- and 

§ Note that many authors call this singularity "central" , especially when studying the collapse 
of massive objects, regardless of its causal characterization. However, this terminology is hard to 
justify in general. For instance, the big-bang singularity in a Robertson- Walker model is simply a 
particular case of a R = singularity, but it would hardly be called central. 
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Vaidya's solution when there is a linear mass function [12J [13J ) , while in other cases 
there is not analytical compactification and the alternative method of studying the 
radial null geodesies has been used in order to get the local causal characterization 
(for example, Vaidya's solution in the general imploding case |14l[|[| or so many 
different collapsing stellar models, like those found in [15J[16J[17J). 

In addition to the analysis of particular cases, some general approaches for 
studying zero-areal-radius singularities have also been carried out by analyzing the 
properties of the radial null geodesies. In particular, it has been shown [TH] that a 
R = singularity is spacelike (and trapped) at a point p in R = if m\ p > and 
timelike (and untrapped) if m\ p < 0, where m = R(l—g' J ' l/ d fl Rd u R)/2. Nevertheless, 
the case m\ p = 0, in which there could be either a regular center or a (spacelike, 
lightlike or timelike) singularity at p, must be analyzed in detail for every particular 
case. The trouble is that this is precisely the most interesting case in many different 
physical situations. For example, what is the outcome of the evaporation of a black 
hole if it gets rid of all of its mass? What transformations can the BH's singularity 
undergo in the process? And, since only particular cases have been treated in 
the current literature [i3j[19j[20j[21j|22j, have all the possibilities been considered? 
On the other hand, if a singularity-free star collapses, can it generate a massless 
singularity at the evaporation event due to the focusing of the different shells that 
constitute its interior? And, in the affirmative case, what can its causal character 
be? 

In order to clarify the importance of these questions let us remind that the 

singularity theorems [5] [UJ show that, given a few reasonable assumptions, a collapse 

can terminate in a gravitational singularity. However, the theorems do not inform 

us about many properties of the singularities |5J. Among others, there is lack of 

information about the type of singularity, the divergence of the energy density 

of matter fields and whether the singularity is hidden from outside view by the 

|| Incidentally, the reader can analyze this case, in which a singularity develops from a regular 
R = with m = 0, to verify that the approach of studying the causal character of n = dR (or 
" g RR ", provided the metric is given in suitable coordinates) in order to characterize a R = 0- 
singularity on the boundary of the spacetime (!) is not a reliable method. 
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formation of a black hole. With regard to the last point, we must emphasize that the 
theorems do allow for the possibility that the singularity could be seen by observers 
close enough to the singularity (in which case we say that there is a locally naked 
singularity) or from the future null infinity (globally naked singularity). In any 
case, there would be a hypersurface (the Cauchy Horizon) beyond which General 
Relativity loses its predictability. The question on whether General Relativity 
contains a built-in safety feature that precludes the formation of naked singularities 
in generic gravitational collapses was put forward by Penrose in 1969 [23] and gave 
rise to what is known as the cosmic censorship conjecture (CCC). Clearly, the causal 
character of the singularity is central in the resolution of this conjecture since, by 
construction, there are always null geodesies which are past incomplete whenever 
the spacetime possesses timelike or past null singularities (see, for example, section 
[2j. Therefore, timelike and past null singularities are always naked. 

Some counterexamples to the cosmic censorship conjecture have been proposed. 
An outstanding collapsing and radiating model can be found in the work by 
Demianski and Lasota [23]. Even if it was not first proposed to be such a 
counterexample, but as an evaporating model, it was later shown [22] that it 
possesses an instantaneous naked scalar curvature singularity at the evaporating 
event. Another counterexample of historical importance was discussed in [T2][T3j[26j 
considering the collapse of null dust modelled by using Vaidya's solution. It was 
shown that it suffices that pure radiation (or null dust) with a sufficiently weak 
wave travelling into an initially flat space-time focuses in R = in order to create a 
null singularity which is at least locally visible. We will not intend now to exhaust 
all the different counterexamples to the cosmic censorship conjecture that can be 
found in the literature, but just to mention some of them we point out that naked 
singularities in spherically symmetric models are also possible in the collapse of dust 
[13] [IB] [27] [28] . perfect fluids [29], general fluids [30], massless scalar fields [31] and 
even in higher dimensions [32] • On the other hand, more general studies on the 
formation of naked singularities in spherically symmetric spacetimes based on the 
study of the radial null geodesies can be found in [33], [34] and, by using ad hoc 
devised procedure, in [3*5] . 
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Notwithstanding the above (incomplete) list of proposed counterexamples to 
the CCC, the subject is still open. This is so because, on the one hand, any specific 
example is unlikely to be considered generic in some appropriate sense and specific 
examples satisfying the CCC exist for the different matter fields above (see [36] [37] 
and references therein). Moreover, in [3H] the weak version of the CCC [3H] has been 
shown to hold for a wide variety of spherically symmetric coupled Einstein-matter 
systems. On the other hand, some examples possessing naked singularities have been 
shown to have some kind of instability. However, it is in general unclear whether 
the instabilities can really hide the singularity (see [36j[40j[41j[42j and references 
therein) . 

Our aim in this paper is to apply the techniques of the qualitative behaviour 
of dynamic systems to the radial null geodesies around every R = singular point 
in order to study the causal characterization of the R = singularities. This will 
allow us to ascertain the relevant quantities (as well as their associated values) 
that determine whether a singularity is essentially spacelike, lightlike or timelike 
(even if the singularity is piecewise timelike or piecewise spacelike, for the case of 
piecewise C 1 boundaries). We will also find out and explicitly state the limits for 
the applicability of our results coming from our specific approach. We would like 
to remark that ours is a geometrical approach requiring only the existence of a 
spacetime, but not the fulfillment of Einstein's equations. Thus, we just try to 
discover the possibilities allowed by this geometrical approach which includes the 
classical as well as the semiclassical framework. With the obtained information 
we will be able to study different possibilities for the final outcome of Black Hole 
evaporation and the different options for the generation of naked singularities. We 
will emphasize new models and possibilities that have not been taken into account 
so far. 

The paper has been divided as follows: In section [2] we study the relationship 
between the null geodesies in the physical and the unphysical spacetimes and we 
revise how to extract information about the causal characterization from them. In 
section [3] we establish a general spherically symmetric spacetime and the conditions 
required for it to have a R = scalar curvature singularity as well as the equations 
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governing its radial null geodesies. Section [4] is an application for the trivial and 
well-known case of a non-zero function m as R tends to zero. In sections [5] and [6] we 
use the theory of qualitative behaviour of dynamic systems to deal with the analysis 
of isolated m(R — > 0) = points, since, as we will see, they turn out to be the 
isolated critical points of the system of differential equations describing the radial 
null geodesies. Specifically, the hyperbolic and the non-hyperbolic cases are treated 
in sections [5] and [6j respectively. In order to exhaust all the different possibilities, 
the non-isolated m(R — >• 0) = points are treated in section [7j The last section 
is devoted to the consequences and applications of our results. In particular, the 
general cases of Black Hole evaporation and the generation of naked singularities 
are analyzed. 

2. Null geodesies and causal characterization 

Let us assume that we are given an oriented two-dimensional Lorentzian manifold 
and that we embed it into an unphysical two-dimensional Minkowskian spacetime, 
as explained in the introduction. Provided that the singular boundary in the 
Minkowskian spacetime is C 1 at a point p we will be able to compute the tangent 
vector to the singular boundary at p. Then, by definition, the causal character of 
the singular boundary at p coincides with the causal character (spacelike, timelike or 
lightlike) of its tangent vector. Furthermore, inspired by the definitions appearing 
in |33], we will also specify that there is a past spacelike (or past lightlike) singularity 
at p if only past-directed causal curves end up at a spacelike (or lightlike) singularity 
at p (see, for example, figures [T](i) and[Tjlv), respectively). Likewise, we say that 
there is a future spacelike (or future lightlike) singularity at p if only future-directed 
causal curves end up at a spacelike (or lightlike) singularity at p (see, for example, 
figures [ljn) and[T|v), respectively). 

Let us now consider a sufficiently small simply connected open set Um in the 
two-dimensional Minkowskian spacetime with corresponding points in the physical 
spacetime (specifically, in an open set IA — defined in section [3] — ) and such that a 
part of its boundary consists of an open interval of the singular boundary containing 
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p. Two nonvanishing null vector fields may be defined on Um such that they are 
linearly independent at each point of Um [E] • The integral curves of the two null 
vector fields provide us with two families (J-\ and J- 2) of (non- necessarily affine 
parametrized) null geodesies. 

2.1. Relating the causal character of singular points to null geodesies 

If we parametrize the C 1 curve describing the singular boundary with a parameter 
A, then the value of the norm of its tangent vector will be given by a continuous 
function A/"(A). The continuity of A/"(A) implies that, if it is bigger than zero at a 
point p corresponding to a certain A = Ao, then it will be bigger than zero for an 
open interval A_ < Ao < A+. So that there is an open interval around p on the 
singular boundary where the boundary is necessarily spacelike. On the other hand, 
in this situation there are only two possibilities for the behaviour of the light-like 
geodesies at p that traverse Um- 

• There is a past ( or future ) spacelike singularity at p =^ A null geodesic of every 
family leaves (resp. reaches) p. (See fig. [l]-i (resp. fig. [l}ii), where we have 
drawn the singularity horizontall}|5} 

Similarly, if A/"(Ao) < then M will be less than zero for an open interval 
A_ < A < A + and there is only a possibility for the behaviour of the light-like 
geodesies at p that traverse Um- 

• There is a timelike singularity at p =>• A null geodesic of one family reaches p 
and a null geodesic of the other family leaves it. (See fig. [l}iii). 

It is useful to note that if one, and only one, null geodesic leaves or reaches 
a point p in the boundary then, by process of elimination, the singularity must be 
lightlike at p. Specifically, there are two cases: 

- If only a null geodesic leaves (or reaches) p then there is a past (resp. future) 
lightlike singularity at p. (An exampld^of this situation is shown in fig. [T]-iv 

% The situation does not change if we draw the singularity with an inclination bigger than -45° 

but less than 45° around p, as is guaranteed by Af(A- < Xq < A+) > 0. 

+ Note that, in general, the singularity does not have to be lightlike all around p. 
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2.2. Intervals of piecewise constant causal characterization 

Let us now classify the intervals in the singular boundary according to the behaviour 
of the null geodesies. (We enumerate them according to the number of geodesies 
that leaves or reaches every of their points): 

1~ ( + )) Intervals where, for every of their points, only a null geodesic leaves (resp. 

reaches) it. Then, according to our previous subsection, the interval is a past 
(resp. future) lightlike singularity. 

2°) Intervals where, for every of their points, a null geodesic of one family reaches 
the point and a null geodesic of the other family leaves it. Then we have seen 
in the previous subsection that the interval can contain points where it would 
be timelike (Af(\o) < 0) and, as a result, subintervals around these points 
would be timelike. The possibility of having a point in this singular interval 
where the singularity is spacelike is trivially discarded since the behaviour of 
the null geodesies does not correspond with the expected one according to the 
first two items in the previous subsection. The possibility of having a lightlike 
subinterval is also discarded since, if we could draw the subinterval, it would 
have an inclination of ±45°, so that its points could not be reached by null 
geodesies from both families as demanded in this case. However, there is still 
the possibility for the interval to have lightlike points only if every neighborhood 
of the lightlike point contains timelike subintervals (an example is shown in 
figure [l}vi). In this case we will say that the interval is piecewise timelike^ 

2~ ( + )) Intervals where, for every of their points, a null geodesic of every family leaves 

(resp. reaches) it. Then we have seen in the previous subsection that the 

interval can contain points where it would be spacelike (Af(\o) > 0) and, as a 

result, subintervals around these points would constitute a past (resp. future) 

* In fact, it can be argued that the distinction between strictly timelike and piecewise timelike 
intervals is not very relevant from a physical point of view. For example, if we are worried about 
the formation of a Naked Singularity in a particular model, we will have one in both cases. 
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Figure 1. The behaviour of null geodesies (arrows) at p in an open C 1 interval 
of the singular boundary (wavy line) is shown for (i) a past spacelike singularity 
at p, (ii) a future spacelike singularity at p and (iii) a timelike singularity at p. 
When only a null geodesic leaves or reaches p there are two possibilities: (iv) p is 
part of a past lightlikc singularity or (v) p is part of a future lightlike singularity. 
In (vi) there is a lightlike singularity at p while the rest of the interval is timelike. 
Thus, it is an example of a piecewise timelike interval. In (vii) there is a lightlike 
singularity at p while the rest of the interval is spacelike. This is an example of a 
piecewise past spacelike interval. 

spacelike singularity. Following the reasoning in the previous item, there is 
only the additional possibility of having lightlike points if every neighborhood 
of the lightlike point in the interval contains spacelike subintervals (an example 
is shown in figure [l}vii (resp. its time reversal)). In this case we will say that 
the interval is a piecewise past (resp. future) spacelike singularity. 

Here finishes our list of possible C l intervals according to the behaviour of 
the two families of null geodesies traversing Um- (Note that there cannot be C 1 
intervals such that every point is reached or left by a total of three or four null 
geodesies traversing Um)- 

As some exact solutions found in the literature show, we could also deal with 
just piecewise C 1 boundaries. In these cases, there would be points where it will 
not be possible to define a tangent vector to the singular boundary and, strictly 
speaking, no matter the chosen approach to ascertain the causal character of the 
boundary we should resign ourselves to deal with piecewise causal characterizations. 
An additional difficulty appears when one analyzes the null geodesies for these 
piecewise C 1 boundaries in order to get its causal characterization. Consider, for 
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example, an open interval in the singular boundary where two null geodesies leave 
the boundary at every of its points. Then we cannot tell with just this scrutiny 
of null geodesies whether there is some point on the interval where the singular 
boundary is not C 1 (see figure [2]). Likewise, a point in the singular boundary 
where the boundary is not C 1 can pass unnoticed either in intervals where two null 
geodesies leave the boundary at every of its points (imagine the time reversal of 
figure [2| or in intervals that are reached and left by one null geodesic of every family 
on every point. However, intervals where, at every point, only a null geodesic leaves 
it (past lightlike singularities) or reaches it (future lightlike singularities) are clearly 
C l (they are defined by ±45°-inclined straight lines) and their causal character is 
therefore strictly determined. 

Let us now call transition points those points in the singular boundary where 
two different intervals from the ones defined above join (for example, a 2° interval 
is followed by a 2~ interval in the singular boundary). This definition implies 
either that the singular boundary is lightlike at a transition point p or that no 
tangent vector to the boundary can be defined at p. Since every piecewise C 1 
boundary will be composed of intervals with a (piecewise) single causal character 
joined through transition points, we will be able to sketch the causal character 
of the singular boundary by identifying the causal character of every (piecewise)- 
single-causal-character interval and overlooking whether the transition points really 
are lightlike. In order to exemplify this, in figure [3]) we have sketched a spacelike 
open interval joined to a timelike open interval through a transition point q. In 
this spirit, from now on we will omit the adverb piecewise in single-causal-character 
intervals. 

2.3. Translation into a physical spherically symmetric spacetime 

Let us assume now that we are working with a spherically symmetric spacetime and 
that we have identified the physical two-dimensional lorentzian surface orthogonal 
to the 2-spheres. Assuming that the invariant areal radius R can be used as 
a coordinate in a local chart of the lorentzian surface, we will work in local 
coordinates {R,u}, where u is a lightlike coordinate. Note that u =constant has a 
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Figure 2. We show, to the left, a proper spacelike open interval of the singular 
boundary and, to the right, a piecewise spacelike singularity (since the -angular- 
boundary it not C 1 at p). Both possibilities are indistinguishable if we only take 
into account that a null geodesic of every family leaves every point in the boundary. 



Sketch 




Figure 3. In a) we draw an open interval of a singular boundary. In b) we show a 
sketch of its causal character. With a scrutiny of null geodesies we can infer that 
(p, q) is a (piecewise) past spacelike interval and (q, r) is a (piecewise) timelike 
interval. However, even if the boundary is lightlike at point q, we can not infer 
this with a simple scrutiny of null geodesies, since the same re-count (one reaches 
q and two leave it) could be obtained if the boundary were C° at q. Therefore, in 
our sketch wc overlook the possible lightlike character of the transition point q. 

clear geometrical meaning, even if no further constraint is imposed on u, since it 
corresponds to an invariantly defined radial null geodesic. Moreover, according to 
our project, with this choice of coordinates a family of radial null geodesies is given 
(u =constant) and we only need to compute the other family of radial null geodesies 
in order to obtain the sketched causal characterization of any R = singularity. 
This is so because, as explained in the introduction, we know that for every null 
geodesic in this physical two-dimensional spacetime there is a corresponding null 
geodesic in the unphysical two-dimensional minkowskian spacetime and viceversa. 
The local translation from the physical (U) to the unphysical (Um) spacetime is 
rather straightforward. We have listed the possibilities for single-causal-character 
open intervals in figure [4j The only noteworthy case appears when an infinite number 
of null geodesies from the second family tends to (or departs from) a singularity at 
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Figure 4. Translation table between the chart in the physical spacetime and the 
sketches in the unphysical spacetime for single causal character open intervals. In 
cases (i) to (iii) at every point in a whole open interval of it's with R = a null 
geodesic of every family tends to (or departs from) the point. (We depict only 
those tending to (or departing from) a point p {R = 0,u = uq} in the interval). 
In cases (iv) and (v) only a family of null geodesies tend to (or depart from) the 
points in a whole open interval of it's with R = 0. In the last possibilities, (vi) and 
(vii), an infinite number of null geodesies (in the grey zone of the R — u diagram) 
tend to a single point {R = 0,u = Uq}, what must be interpreted as a lightlike 
singularity in the unphysical spacetime. 

{R — 0, u — Mo}, which corresponds to the existence of a lightlike singularity, as can 
be seen in figures [3]- (vi) and |3]-(vii). 

3. General spherically symmetric metric 

Let us consider an oriented four-dimensional spherically symmetric space-time V. 
In order to study the local causal behaviour around R = we will consider a local 
chart endowed with coordinates = {u, R, 9, <p} (« = 0, 1, 2, 3) and an open set 
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U = {(u,R)\ — 5u < u < Su, < R < 5R}. The line-element can be expressed in 
this local chart as 

ds 2 = -e^xdu 2 + 2ee w du dR + R 2 dtf, (1) 

where R is the areal radius, e 2 = 1, j3 and x are assumed to be at least C 2 ~ function^ 
on {u, R} in the local chart, (3 is also assumed to be bounded as it approaches R = 0, 
and dQ 2 = d6 2 + sin 2 6d<f 2 . 

In spherical symmetry one can define the scalar invariant 

m =^(l-g^d,Rd u R) 

(see [45J and also [TS] , [IB] and references therein) . If one computes m for the metric 
([I]) it is easily checked that 

m = |(l- X ). (2) 

On the other hand, if we do not want R = scalar curvature singularities, the 
scalar invariants polynomial in the Riemann tensor must remain finite at R = 0. 
It is well known that there are only four algebraically independent scalar invariants 
associated with a general spherically symmetric metric [J7]. We can take, for 
example, pE] pE] [50] : 

K = g aj gP 5 R afhS , 

where S a ^ = Rj 3 — d^TZ/ '4, being R a ^ the Ricci tensor and 1Z the curvature scalar; 

Ca^s the Weyl tensor, C a p^s = (C a p 7 s+i *C a ^s)/2 is the complex conjugate of the 

selfdual Weyl tensor being *C a/ 3 7( 5 = € a/3fJjU C flu ^ s /2 the dual of the Weyl tensor. If one 

evaluates these invariants for Q using ^ one arrives to the following statement 

[50j[51j: All scalar invariants polynomial in the Riemann tensor will be finite at 

jj This is a minimum requirement which guarantees the existence and uniqueness of null geodesies 
in U llj. It also enables to use Einstein's equations (in case one works in the framework of GR). 
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R = 0, preventing the existence of scalar curvature singularities if, and only if, 

lim m = lim — = lim — - = lim — — — = 0, 

fl->0 r-^o R R->0 R 2 b^o R 

lim = &(«), lim^ = m3 H (3) 

where (3 (u) = Km^o /?(«, R) and both (3 2 (u) and m 3 (u) are finite functions. 

Moreover, if the set of scalar invariants is finite at R = there will be only two 
algebraically independent scalar invariants at R = 0, say 1Z and n, since W2 = 
and 343rx 6 = 3087r 2 4 . 

Let us finally remark that one of the interesting properties of the scalar 
curvature singularities is that their existence can not be an artifact of the coordinate 
system used. We now know that if conditions ^ are not fulfilled then for any curve 
approaching R = there will be at least one scalar invariant that will grow without 
limit along the curve as it approaches the R = 0-singularity. It is obvious that, if 
one uses different local coordinates, the value of the scalar invariant at every point 
of the curve must coincide, when computed with the new local coordinates, with the 
values obtained at the same points with the coordinate system used in this article. 
Therefore, no matter what local coordinates are used, the same scalar invariant will 
grow without limit along the curve, what unmistakably identifies a R — 0-scalar 
curvature singularity. 

3.1. Radial null geodesies 

We choose u growing to the future. Then 

1=- = --— (4) 
dt 2RdR K 1 

is a future directed radial null vector tangent to the null geodesies u =constant 
(namely, the T\ family) and i is a future directed parameter. Since the expansion 
[TTj of these null geodesies is given by 

if £ = —1 (or +1), the expansion is positive (negative, respectively) and, they are 
outgoing (ingoing, respectively) radial null geodesies directed towards increasing i?'s 
(decreasing i?'s, respectively) according to Q. Therefore, the behaviour of the T\ 
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family, with regard to whether the null geodesies are coming or are directed towards 
the R = 0- singularity, is absolutely defined by the sign of e, a fact that we will use 
throughout this article. On the other hand, 

is a future directed radial null vector such that 1 • k = —1 and k is a future directed 
parameter. For later purposes let us write explicitly the equations governing the 
geodesies that have k as its tangent vector field (the Ti family): 

e[R-2m(u,R)} 

(6) 

In this case the expansion is given by 
e k = 2 ^(R-2m). 

If at a given 2-sphere x > R > 2m) and e = —1 (or +1), k is tangent 
to a family of null geodesies with negative (positive, respectively) expansion, the 
areal coordinate R decreases (increases, respectively) along them according to ^ 
and, therefore, these radial null geodesies are ingoing (outgoing, respectively) in the 
considered 2-sphere. However, it is interesting to note that if, in a given 2-sphere, 
X < and e = —1 (e = +1) then the two radial null vectors have both positive 
(negative) expansion which means that the 2-sphere is a closed surface trapped to 
its past (future, respectively). 

4. Characterization of the singularity: Case lim,^^ m(u, R) ^ 

Let us consider, in U, the set of radial null geodesies belonging to the family Ti that 

approach R = (either toward their past or their future). For every such geodesic 

we define k to be the value of its parameter such that lim re _>. reo (i2(«;), u(k)) = (0, uo), 

where the precise value uq G {—Su, 8u) depends on every geodesic. Now, if 

lim K ^. K0 m(u, R) ^ along every of these null geodesies, according to rt6J), 
dR 

lim — — = — 2e lim m(u, R) (7) 

K— >H (IK, ft— >K0 

du 

hm — = 

k->«o du 

What implies that: 



Causal characterization of singularities 17 
• If lim K _> Ko m(u, R) > in U then the two different families of radial null 



geodesies are outgoing (e = —1) or ingoing (e = +1 ) in U (see sect. 3.1 for 
the u =constant family) and, therefore, the singularity is spacelike at R = 0. 

• If \im K ^ Ko m(u, R) < in U then one radial null geodesic is outgoing while the 
other is ingoing and, therefore, the singularity is timelike at R = 0. 

As a corollary we have the following more workable and well-known [IS] result: 
if lim^j^Q) m(u, R) exits and it is positive for all uo G U then the singularity is 
spacelike at R — and if lim( u _>.o,.R->o) m(u, R) exists and it is negative for all u G W 
then the singularity is timelike at R = 0. 



5. Characterization of the singularity: Case m(u = 0, R = 0) = and 

rh, u {u = 0,R = 0) ^ 

In this section we will assume, as usual, that m and /3 are C 2 ~ in U and also that C 1 
auxiliary extensions (devoid of any physical meaning) fh and /3, respectively, exist 
in U = {(«, R) \ — 5u < u < Su, —5R < R < d~R} such that, formally, fh = m and 
P = P for {-5u <u <5u,0 < R< StiQ 

In the case considered in this section (m(u = 0, R = 0) = and m, u (w = 0, i? = 
0) 7^ 0) there is a scalar curvature singularity at R = u = since the regularity 
conditions ^ are violated. 

Without loss of generality, we will consider from now on j3(u, 0) = since, if 
it was not, we can always perform a coordinate change u — > v! such that the new 
coordinate u' were defined by du' = e 2/3( > u ' * > du. 

Now u = R = is a critical point of system ([6]). In order to analyze the 
qualitative behaviour of the radial null geodesies we try a linear approximation at 
the critical point: If we define f R = e[R — 2m(u, R)] and f u = 2Re~ 2/S , then the 
linearization matrix is = f a ,p (u — 0, R — 0): 

I e[l-2m, B (0,0)] -2em, tt (0,0) 
I 2 

f Note that in the literature one usually finds the same functions m and /3 used for R < 0, what 
could be considered as the natural extension. 
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The characteristic roots for this matrix are: 



s(l - 2m, R (0, 0)) ± ^(1 - 2fh, R (0, 0)) 2 - 16em, u (0, 0) 



(9) 



The critical point u = R = will be hyperbolic [6J if none of the characteristics 
roots have zero real part. In order to analyze this for our case, let us define 
A = (1 - 2m, R (0, 0)) 2 - 16em, u (0, 0). Then u = R = will be hyperbolic if 



and, in case A < 0, if the extra-condition m,R (0, 0) 7^ 1/2 is satisfied. Under the 
assumption that m and (3 are C 1 we can apply the Hartman-Grobman Theorem 
[S]. According to this theorem, if the critical point is hyperbolic then the behaviour 
of the nonlinear system near the critical point is qualitatively determined by the 
behaviour of the linear system x a = A^x 13 . Furthermore, the theorem implies that 
the qualitative behaviour is independent of the extension chosen for m and /3 when 
{—5u < u < 5u, —5R < R < 0} as long as it is a C 1 extension. In other words, and 
taking into account (j8J), the results will only depend on the partial derivatives of m 
as we approach R = u = 0. Note that in the hyperbolic case and near the critical 
point the curves defined by f R (u,R) = and f u (u,R) = behave approximately 
as straight lines in the {u, R} plane intersecting at the point u = R = 0, which is 
an isolated critical point 0. 

5.1. Different possibilities 

• If A < and m,R (0, 0) ^ 1/2 the characteristic roots are complex conjugate of 
each other with nonzero real part. This means that the critical point u = R = 
is a focus |6j|8j. In this way, the future directed^ radial null geodesies in a 
neighborhood of the critical point must start at R = 0, u < and end at 
R = 0, u > after going around the critical point. Therefore, 

— If e — — 1 then the R = singularity must be spacelike near the critical 

point for u < while it must be timelike near the critical point for u > 0. 

I Throughout the article we will use the fact that the radial null geodesies are future directed 
in order to interpret in correct physical terms the results provided by the qualitative theory of 
dynamic systems. 



m m (0,0) ^0 



(10) 
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Figure 5. In the left of this figure we sketch the center-focus behaviour of the 
second family of radial null geodesies in the e = — 1 case near u = R = when 
det^4 > and A < 0. (In fact, it is only a half of the center-focus behaviour 
since we are only interested in R > 0). Note that the sketchy representation gives 
rise to the illusory discontinuousness of some curves. In the right we show the 
corresponding sketched Penrose diagram around u = R = with its characteristic 
spacelike-timelike behaviour in the s = — 1 case. 



— If e = +1 then the singularity must be timelike near the critical point for 
u < while it must be spacelike near the critical point for u > 0. 

• If A < and m,#(0,0) = 1/2 the point is not hyperbolic. However, even if 
the linearization does not suffice to distinguish the exact qualitative behaviour 
it can be guaranteed that u = R = will be either a focus or a center jH] [7]. 
Therefore, the same characterization than in the above item applies. 

• If A > and erh, u (0, 0) > then the characteristic roots are real, with the 
same sign and distinct so that the critical point is a node [6 J \S\. The slope of the 
radial null geodesies ending or starting at u = R = is £ = h m K-s-K u(k)/R(k). 
Using g we find £ = 2e/(l - 2{m, R (0, 0) + £ rh, u (0, 0)}). The two real roots 
£o± of this quadratic equation are 



— If T = e(l — 2m, R (0, 0)) > then \± > 0, £ ± > and the node is unstable. 



See figure [5] 



(1 - 2m, R (0, 0)) ± y/(l - 2m, R (0, 0)) 2 - lQem, u (0, 0) 

4rh, u (0,0) 



.(11) 
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Figure 6. In the left of this figure we sketch the unstable node behaviour of the 
second family of radial null geodesies in the e = — 1 case near u = R = when 
dct^4 > 0, T > and A > 0. The two straight lines crossing at u = R = indicate 
the critical directions of the radial null geodesies starting at u = R = 0. From 
these geodesies only one leaves the critical point with direction £o+ while all the 
others leave it with direction £o- • In the right we show the corresponding Penrose 
diagram around R = u = with its characteristic R = u = O-lightlike singularity. 
Here we have pointed out the single radial null geodesic that leaves the critical 
point in the direction £o+- 

Then it can be shown [7] that around u = R = all, except for one, radial 
null geodesies must start at this point with a definite slope which is £ + 
for e = +1 and £ - for e = —1. The exception is just one null geodesic 
starting with slope £o- f° r £ = +1 and £o+ for £ = —I- The behaviour of 
the second family of null geodesies (J-jj) for the e — — 1 case together with 
the corresponding sketched Penrose's diagram is shown in figure [6j 
— If T < then A± < 0, £o± < and the node is stable. The comments for 
this case are similar to the case above. We summarize the behaviour of 
the second family of null geodesies for the e — — 1 case together with the 
corresponding sketched Penrose's diagram in figure [7) 

• If A = and em, u (0, 0) > then the characteristic roots are real and equal (A) 
so that the critical point is a degenerate node |6j[8j. There also exists an only 
critical direction (£ )- 
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Figure 7. In the left of this figure we sketch the stable node behaviour of the 
second family of radial null geodesies in the e = — 1 case near u = R = when 
detA > 0, T < and A > 0. In the right we show the corresponding Penrose 
diagram around R = u = with its characteristic R = u = O-lightlike singularity. 



— If T > then A > 0, £o > and the node is unstable. Then it can be 
shown [7] that around u = R = all radial null geodesies must start at this 
point with a definite slope which is £o f° r e — ±1- Therefore, the result is 
similar to figure [6] with equal critical directions. 

— If T < then A < 0, £ < and the node is stable. Therefore, the result is 
similar to figure [7] with equal critical directions. 

• If A > and em, u (0,0) < (no matter the value of e(l — 2m, #(0,0))) 
the characteristics roots are real with opposite sign and the critical point is 
a saddle [E][B]- Likewise, the two critical directions also have opposite sign 
since sign(£o-) — e — — sign(£o+)- The behaviour of the second family of radial 
null geodesies for the e — — 1 case together with the corresponding sketched 
Penrose's diagram is shown in figure [8} 

We have collected the results for rh, u (u = 0, R = 0) ^ in figure [9j 
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Figure 8. In the left of this figure we sketch the saddle behaviour of the second 
family of radial null geodesies in the e = — 1 case near u = R = when detA < 
and A > 0. In the right we show the corresponding sketched Penrose diagram 
around u = R = with its characteristic timelike-spacelike for the e = — 1 case. 




Figure 9. Characterization of the singularity when m(u = 0, R = 0) = and 

m, u (u = 0,R = 0) ^ 0. 
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6. Characterization of the singularity: Case m{u — 0, R — 0) = and 

rh, u {u = 0,R = 0) = 



In this case u = R = is an isolated critical point of system (|6j), but, as we showed in 
the previous section, since m, u (u = 0, R = 0) = the critical point is not hyperbolic. 
In fact, if we demand fh,R (0, 0) ^ 1/2 the critical point will be semi-hyperbolic 
and we will have to use a different approach. (The theory regarding the qualitative 
behaviour of these points can be found in [U] PH] • The reader can also find a summary 
of the main results in appendix A). In order to apply the theory for semi-hyperbolic 
points we must now demand, first, the existence of a natural number n > 2 such 
that lim^o ^o) d n m/du n (u, R) ^ 0, while lim^^R-^Q) d % m/du l {u, R) = for 
2 = 1, ...,n — 1, and, second, of C n extensions^] in and {3. 
First we define 

2e 

x= ; -R + u (12) 

1-277^(0,0) 1 ; 

y = R (13) 
t = e{l - 2m, R (0, 0))«. (14) 
Then the system ^ can be rewritten in the normal form 

f = P2{x,y) 

where P 2 (x, y) = P 2 (u(x, y), R(y)), Q 2 (x, y) = Q 2 (u(x, y), R(y)), 



P 2 (u, R) = 2 % , ( Re-^ _ R-*m{u ; R) 

V ' l-2m,«(0,0) V l-2m JjB (0,0. 

Q 2 {u,R) = 2 



0) 

m, R (0, 0)R — m(u, R) 



l-2m, fl (0,0) 

and they satisfy P 2 (0,0) = Q 2 (0,0) = P 2 , x (0, 0) = Q 2 , x {0,0) = P 2 , y {0,0) = 

Q2, y (0,0) = 0. The assumed degree of differentiability together with the implicit 

function theorem guarantees that the equation y + Q 2 (x,y) = has a solution 

y = (p(x) which can be written (if one considers the Taylor polynomial for fh(u, R) 

§ Private communication with the authors of [TU], where C°° is assumed independently of the 
value of the finite n. 
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and does some algebra) as the truncated series expansion[jj] 
2 0^/0u"(O,O) 

y ~ n\l-2m, R (0,0f + -' (15) 
where n > 2 is assumed to be finite. If we define the function ip{x) = P2(x,(p(x)) 
then its truncated series expansion will have the form 

i>{x) = A n x n + (16) 

where, in our case, 

f) n ffi 

A„ = ec 2 — (0,0), (17) 

and c 2 = 4/(n!(l - 2m, R (0, 0)) 2 ) > 0. 

The theory of the qualitative behaviour of dynamical systems tell us that we 
only need to know whether n is even or odd and the sign of A n to tell the qualitative 
behaviour of the geodesic curves around the critical point which can now be a saddle 
node, a topological saddle or a topological node. In fact, the saddle node behaviour 
is the only possibility that we have not treated yet. While all saddle nodes are the 
union of one parabolic and two hyperbolic sectors [9], we still have to study the 
geodesies only in U, what implies R > 0. There appear four different possibilities 



that we show in the self-explanatory figures 10 11 12 and 13 for the case e = —1. 
We have collected all the possibilities for m, u (u = 0, R = 0) = and 



m, R (0, 0) ^ 1/2 in figure 14 



Finally, let us analyze the case with m, u (0, 0) = and m,n (u — 0, R — 0) = 
1/2. In this case the critical point will be nilpotent [10]. Again we refer the reader 
to [H] PH] for the theory of the qualitative behaviour of these points and to appendix 
B for a summary. In order to apply the theory of nilpotent points we must demand 
the existence of, first, a natural number n > 2 defined by n = Max{i,k}, where 
i is the lowest value of j such that lini( M _>o,_R^o) d^m/du^u, R) ^ and k = I + 1, 

|| Note that, the requirement on the degree of differentiability for m and (3 -taking into account 
Taylor's theorem- allows us to write some functions derived from m and j3 as a Taylor polynomial 
plus a remainder term. From now on we will only write the first non-zero term of the Taylor 
polynomial (the only that is guaranteed to exist thanks to our assumptions) and suspension points 
-as in ( 15 ) . 
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Figure 10. In the left of this figure we sketch the behaviour of the second family 
of radial null geodesies in the e — — 1 case near u = R = when fh, u (0,0) = 0, 
A„ < with n even and fh,B. (0,0) > 1/2. Only the R > part of the saddle 
node structure is shown. The straight line marks the regular critical direction £ 
[8] for the curve starting at u = R = that separates a hyperbolic sector (H) from 
the parabolic sector (P). In the right we show the corresponding Penrose diagram 
around u = R = with its spacelike singularity. 
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Figure 11. In the left of this figure we sketch the behaviour of the second family 
of radial null geodesies in the e — — 1 case near u = R = when fh, u (0,0) = 0, 
A„ < with n even and m,R (0, 0) < 1/2. We also show by means of a straight line 
the regular critical direction £o- In this case the two hyperbolic sectors (H) and 
the parabolic sector (P) are partially visible for R > 0. In the right we show the 
corresponding Penrose diagram around u = R = with its lightlike singularity. 
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Figure 12. In the left of this figure we sketch the behaviour of the second family 
of radial null geodesies in the e — — 1 case near u = R = when fh, u (0,0) = 0, 
A„ > with n even and fh,R (0, 0) > 1/2. The two hyperbolic sectors (H) and the 
parabolic sector (P) are partially visible for R > 0. We also show by means of a 
straight line the finite critical direction £q- I n the right we show the corresponding 
Penrose diagram around u = R = 0, which is a lightlike singularity. 
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Figure 13. In the left of this figure we sketch the behaviour of the second family 
of radial null geodesies in the e — —1 case near u = R = when fh, u (0,0) = 0, 
A„ > with n even and fh,n (0,0) < 1/2. We also show by means of a straight 
line the regular finite critical direction £o of the curve starting at u = R = that 
separates an hyperbolic sector (H) from the parabolic sector (P) for R > 0. In 
the right we show the corresponding Penrose diagram around u = R = with its 
timelike singularity. 
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Figure 14. Characterization of the singularity when m(u = 0, R = 0) = 0, 
m, u (u = Q,R = Q)=0 and m, R (0, 0) 7^ 1/2. 
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where I is the lowest value such that lim( u ^ ,R^o) d l (dm/dR)/du l ^ and, second, 
C n extensions^] rh and (3. 

In this case we simply define x = u/2, y = R, t = n in order to rewrite the 
system (|6|) in the normal form 

y I I'Ar.f,) 



dt 

i = Q2(x,y), 



where 



P 2 {x,y) = P 2 (u{x),R(y)), P 2 {u,R) = R(e~ 2 ^ - 1), 
Q 2 (x, y) = Q 2 (u(x), R(y)), Q 2 (u, R) = e{R - 2m(u, R)) 

and they satisfy P 2 (0,0) = Q 2 (0,0) = P 2 ,x(0,0) = Q 2 , x (0,0) = P 2 , y (0,0) = 
Q 2 , y (0, 0) = 0. The solution of the equation y + P 2 (x, y) = which, by the implicit 
function theorem in general takes the form y = <p(x), in this case is simply y = 0. 
Following the general procedure, If we define the function ip(x) = Q 2 (x,(p(x)) then 
the assumed degree of differentiability guarantees that its truncated series expansion 
will have the form 

ip{x) = a k x k + ... (18) 

where k > 2 is assumed to be finite and, considering the Taylor polynomial for rh, 
one gets 

2 k+1 d k rh, . . , 

°' = - £ ^a^ (0 - 0) - (19) 

On the other hand, if we define a(x) = P 2 , x (x, <f(x)) + Q 2 , y (x, <p(x)) then either its 
truncated series expansion can be written as 

a{x) = b n x n + ... (20) 

where 

2 n+1 d n fdm\ . . 
&n = -e-r^ l -KB (0,0) (21) 



n\ du n \dR J 

if b n 7^ or, alternatively, o~(x) = in which case h n = Vn. 

% Private communication with the authors of [TU], where C°° is assumed independently of the 
value of the finite n. 
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Figure 15. In the left of this figure we sketch the behaviour of the second family 
of radial null geodesies in the e — — 1 case near u = R = when fh, u (0,0) = 0, 
fh,n (0, 0) = 1/2, a,k < with k odd, 6„ > with n odd and either n < (k— 1)/2 or 
n = (k — l)/2 and A = 6^ + 2(fc + l)ofc > 0. Note that we also show the behaviour 
for negative R. The fact that the hyperbolic sector (H) and the elliptic sector (E) 
that constitute the elliptic region appear for R > and for all u in U induces the 
unusual light-like singularity that we have drawn in the right. Specifically, all the 
radial null geodesies (of the second family) that come from the R — 0-singularity 
for < u < —Su must reach in their future R = u = 0, what implies that there is 
a future lightlike singularity at R = u = (see figure EKdi)). On the other hand, 
all the radial null geodesies that reach (R = 0, > u > 8u) come from R = u = 0. 
what implies that there is also a past lightlike singularity at R = u = (see figure 

ER). 



The theory of the qualitative behaviour of dynamic systems tell us that we 
only need to know whether n is even or odd and the sign of and b n to tell the 
qualitative behaviour of the null curves around the critical point which can now 
be a saddle node, a topological saddle, a topological node, a cusp, a focus-center or 
an elliptic region. In fact, only the equilibrium state with elliptic region behaviour 
provide us with a behaviour for the singularity that we have not encountered so far. 
We show this possibility in figure [15] 

We have collected all the possibilities for m, u (u = 0, R = 0) = and 



fh,n(0, 0) = 1/2 in figures 16 and 17 with supplementary details in their 
corresponding captions. 

The cases analyzed so far cover all the possibilities whenever u = R = is an 




Figure 16. Characterization of the singularity when fh(u = 0, R = 0) = 0, 
(u = 0, i? = 0) =0, TO,/{(0,0) = 1/2 and A: is odd. In the extra- conditions 
we have use p = (fc — l)/2 and A = b\ + 2(k + l)ak- The reader can check (see 
appendix B or [21 HD] ) that the first row corresponds to a topological saddle, the 
second to a focus-center, the third and fourth to a topological node and the fifth 
and sixth to an elliptic region. Note that the topological node and the elliptic 
region have two possibilities depending on the sectors appearing in R > 0, what 
depends on the sign of b n . 
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Figure 17. Characterization of the singularity when rh(u = 0, R = 0) = 0, 
(k = 0, R = 0) =0, fh,R (0,0) = 1/2 and fc is even. In the extra- conditions 
we have use p = fc/2. The reader can check (see appendix B or [TU]) that the 
first and second rows correspond to cusps and the rest to saddle nodes. Note that 
both the cusp and the saddle node have different possibilities depending on the 
sectors appearing in R > 0, what now depends on the sign of (and also of b n 
in the saddle node case). 
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7. Characterization when m(—5u < u < 5u, R = 0) = 

In this case the curves f R {u, R) = e[R - 2m(«, R)} = and f u (u, R) = 2Re~ 2 P = 
must cross in an interval made up of points which are non-isolated critical points for 
the system of differential equations (|6|. Assuming the existence of the C 1 extensions 
fh and {3, the equation for the radial null geodesies ([6j can be easily written on the 
interval, up to first order, as 

C W = f( M ,fl = 0)= g(1 - 2 ^'"' 0)) , (22) 

where we have used the fact that rh(—5u < u < 5u, R = 0) = implies that m, u = 
in the interval. Let us consider the characterization of the singularity for the case 
e = — 1 (the case e — +1 can be easily obtained later on as its time-reversal). Then 
the u=constant radial null geodesies are outgoing while 

• If m,R (u, R = 0) < 1/2 in hi the second family of radial null geodesies satisfy 
C < and they are ingoing. Therefore, the singularity there^is time-like. 

• If fh,R (u, R = 0) > 1/2 in hi the second family of radial null geodesies satisfy 
C > and they are outgoing. Therefore, the singularity there is space-like. 

• If rh,R (u,R — 0) = 1/2 in U then ( = 0. In this way, the trajectories of the 
second family of radial null geodesies tend to be parallel to the R = interval 
the closer they are to the singularity and, thus, the singularity is light-like. 



We have sketched the three situations in figure 18 for e — — 1. (It is not necessary to 



sketch the case e — +1 since it produces the time reversal diagrams for every case) 



8. Concluding remarks and some applications 

Let us assume that the reader finds a specific singular model and wants to check 
the causal character of its R = 0-singularity by means of the results provided in this 
article. In order to clarify the applicability of our approach and to make the reader's 
task easier, we would now like to summarize the assumptions that have been made 
+ In case there is a singularity, since this case includes the regular one 
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Figure 18. We have sketched above the u — R graphics for the three different 
cases when there is a non-isolated critical point and e = — 1. The order is, 
from the left to the right, fh,n(u,R = 0) < 1/2, rh,R (u, R = 0) > 1/2 and 
fh,n(v,,R = 0) = 1/2. Below these graphics we have drawn the corresponding 
sketched Penrose diagrams. 

along the article and the path that the reader should follow. On the one hand, we 
have assumed that 

- We are working with a time-orient able spherically symmetric spacetime 
possessing a local chart endowed with coordinates {u, R, 9, ip} (see section [3j. 
In this way, we assume that the invariant areal radius R can be used as a 
coordinate in our local chart. 

- The metric, written as in ([I]), depends on functions (3 and \ which are at least 
C 2 ~ in the local chart (this is a minimum requirement in order to guarantee 
the existence and uniqueness of null geodesies. Note however that the degree 
of differentiability is usually required to be higher). is also assumed to be 
bounded. 

- The local study is carried out for finite values of u and, specifically, we have 
chosen to work around u = (if necessary, it suffices a simple coordinate change 
u = u+ const ant). We assume that, at least, (future or past directed) radial 
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null geodesies of the T\ family (u =constant) reach R = in an interval around 
u = (— 5u < u < Su). 

- The spacetime has a singular boundary that is (in the unphysical spacetime), 
at least, piecewise G x . 

- In the cases where an extension (m) of the invariant m is required, then we 
have not considered working around non-isolated critical points of ^ that are 
accumulation points of the set m(R = 0, u) = 0. 

On the other hand, the reader should be aware that, for every case, our use 
of the qualitative theory of dynamic systems implies that the complete path could 
only be followed if some extra-requirements on the differentiability of the functions 
m and /3 are satisfied. Let us then summarize the aforesaid path along with every 
differentiability requirements: 

• If lim( u ^ UOi R_ s>0 ) rn(u, R) ^ (Vuo G U, i.e., for an interval of iz's) there is not 
assumption required here^] The limit can be either finite or infinite and the 
characterization depends only on whether it is positive (spacelike singularity) 
or negative (timelike singularity) (final corollary in section Q. 

Note that if lim( u _ >UOv R_ s>0 ) rn{u, R) does not exist then we can still evaluate m 
along the radial null geodesies of the second family: lim K _> K0 m(u, R) (see section 
[1]). If it is not zero, no matter if it is finite or infinite, the characterization 
depends only on whether it is positive (piecewise spacelike singularity) or 
negative (piecewise timelike singularity). On the other hand, if only the 
directional limit exists and it is zero the method is not conclusive. 

• If lim^-^fl-^o) m(u, R) = (Vw G U, i.e., for an interval of it's) then we 
demand the existence of C 1 extensions for m and (3. If lim( M _ s . UO) ^_ !> o) m ,R (u, R) 
is less than 1/2 in U then there is a timelike singularity, if the limit equals 1/2 
in U then there is a lightlike singularity and if it is greater than 1/2 then there 
is a spacelike singularity (section [7]). 

• If lini(. u _ > o ) .R->-o) m(u, R) = and lim^^o.ii-^o) m ,u (u, R) 7^ then we demand 
the existence of a C 1 extension for m and /3 (section [H]). In this case the results 

* Except for the obvious minimum degree of differentiability (C 2 ~) for m and /3 in Li. 
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can be found in figure |9j 

• If lim^o.i^o) rn(u, R) = (and only for u -> ), lim (u _ K) ,.R->o) rn, u (u, R) = 
and liHi( n ^ 0jjR ^ ) (w, R) ^ 1/2 then we demand, first, the existence of 
a natural number n > 2 such that lim( u _>.o,fl->o) d n m/du n (u, R) ^ 0, while 
lim( u _ >0 ,ii->-o) d l m/du l (u, R) = for i = 1, ...,n— 1, and, second, of C n extensions 
m and (3 (sectional). In this case the results can be found in figure 



14 



If lim^o^o) m(u, R) = (and only for u -> 0), lim^o^o) rn, u (u, R) = 
and lini( u _ > o,.R-»-o) rn,R (u, R) = 1/2 then we demand the existence of, first, a 
natural number n > 2 defined by n = Max{i, k}, where i is the lowest value of 
j such that lim( u ^ ,R^o) &mfdv?{u, R) ^ and k — I + 1, where I is the lowest 
value such that lim^o ^o) d\dm/dR)/du l ^ and, second, C n extensions fh 



and (3 (section pi). In this case the results can be found in figures 16 and 17 



Assuming that we use a single local chart and that the required C n extensions 
exist for every case, an interesting question appears: Are there any forbidden 
possibilities for the causal behaviour of the singular boundary? If one inspects 
the figures in the article it is easy to see that, indeed, some possibilities do not 
appear. Specifically, there are four possibilities which do not appear due to the 
degree of differentiability required. We have collected them under the name of "135° 
singularities" (see figure [l9](l-4)) . This is expected from the qualitative theory of 
dynamic systems since the theorems involved only allow for specific combinations 
of sectors. On the other hand, the use of a local chart does no allow us to work 
around transition points that can be reached by two radial null geodesies and left 



by another two (pointed singularities, see figure 19 5 and 6)) since, as explained in 



subsection |3.1[ the u =constant radial null geodesies must always be outgoing or 
ingoing in a single local chart. 



8.1. The role of the scalar invariant m in the causal characterization 

As we have seen, despite the function is involved in the behaviour of the radial 
null geodesies ([6]), the causal characterization of the zero-areal-radius singularity 
around u = R = for all the different possibilities depends in quantities (A, T, 
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Figure 19. The forbidden possibilities if one works with a single local chart and 
assuming the existence of the required C n extension. From 1) to 4) they are 
particular transitions from timelike to lightlike singularities and from spacelike to 
lightlike singularities forbidden by the requirement of the degree of differentiability 
for the extension. Specifically, the forbidden possibilities are those that in our 
sketches are represented forming an angle of 135 degrees in the spacetime. From 
5) to 6) we show the pointed singularities that can not be described around p with 
a single local chart in our coordinates. We explicit a possibility showing that the 
u =constant geodesies cannot be always ingoing or outgoing as required (in both 
cases Mi is ingoing while it 2 is outgoing). 

detA, A n , a,}., b n ) that are just obtained from the knowledge of the lowest-order 
non-zero derivatives of m(u,R) as u and R tend to zero. Then we can state the 
following result: 

Proposition 8.1 Provided the listed assumptions are satisfied the knowledge of the 
scalar invariant m is all that is required in order to obtain the sketched causal 
characterization of the R = singularity in a spherically symmetric space-time. 

We can refine this result taking into account that the knowledge of m is not 
a necessary, but a sufficient condition in order to characterize the singularity in 
this approach. By inspecting the quantities defining the character of the singularity 
we see that the knowledge of the lowest non-zero lim u ^o,J?-s>o d k m/du k {k > 0) and 
lim u ^ ,i?^o d n m,ji /du n (0, 0) {n > 0) suffices to characterize the singularity. Then, 
assuming the existence of the extension for m with its corresponding degree of 
differentiability and using Taylor's theorem up to that degree we see that a rough 
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knowledge of fh(u, R = 0) and rh,R (u, R = 0) (note that they are evaluated just at 
R = 0) suffices to characterize any R = singularity in U. 

8.2. Lightlike singularities and shell-focussing nakedness 

When dealing with naked singularities that form due to the gravitational collapse, 
one usually looks for the generation of timelike or lightlike singularities. Provided 
that our assumptions are satisfied, the timelike case is straightforward; it requires 
either lim^o m < in an interval of u's in (section [4]) or lim^om = (if the 
extra condition lim^Qm,^ (u,R) < 1/2 is satisfied for some proper interval of w's in 
U) (section Q. The lightlike case, however, requires a careful study to which we will 
devote this subsection. For example, assuming the existence of the corresponding 
extension for m and by inspecting the hyperbolic and the semi-hyperbolic cases 



(figures M and 14 , respectively) we arrive to the 



Proposition 8.2 Assuming the existence of, first, a natural number n > 

1 such that lim( u _>.o,.R->o) d n m/du n (u, R) ^ (while, in case n ^ 1, 
lim^o^o) d l m/du l (u, R) = fori = l,...,n — 1) and, second, of C n extensions 
rh and (3, if a R = singularity satisfies m(0, 0) = 0, 

e^(0,0)[l-2m, fl (0,0)f +1 > 0, (23) 
[1 - 2m, R (0, 0)f - I6em, u (0, 0) > 0, (24) 

then it has a light-like singularity at R = u = 0. 

Furthermore, if e(l — 2m, # (0, 0)) > there is a past lightlike singularity, 
whereas if e(l — 2m, # (0, 0)) < there is a future lightlike singularity. 

As a particular application let us consider the case of Vaidya's radiating 
solution. Since in this solution m only depends on u we will have 

Corollary 8.1 Any Vaidya metric admitting, first, a natural number n > 1 such 
that lim u ^o d n m(u)/du n ^ (while, in case n ^ 1, lim u ^od l m/du l = for 

2 = 1, ...,n — 1 ) and, second, a C n extension rh, will develop a light-like singularity 

jj Think, for instance, of the timelike singularity in the Reissner-Nordstrom solution, where a 
negative limit is justified due to the presence of electrical charge. 
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at R = u = if 

m{u = 0) = 0, 
dm , , 1 

e *r<» = 0) £ Ye- 

e »^(„ = 0) > 0. 

This corollary generalizes previous results on Vaidya's metric with e = +1 for 
the linear [13] [26] and the non- linear [H] cases. 

Among the lightlike singularities some of them are persistent naked singularities. 
By "persistent" we denote a naked singularity such that a whole family of future- 
directed lightlike radial null geodesies emerges from it. It has been pointed 
out that these naked singularities can appear as a consequence of shell-focusing 
[I3][IS][16j[17j[26j[52] during the collapse and that they are, therefore, relevant to 
the cosmic censorship conjecture. If we want to avoid that the 2-spheres close enough 
to the R = 0-singularity become past-trapped closed surfaces (usually only future- 
trapped 2-surfaces are considered as possible byproducts of collapse) we should also 



demand e = +1 (see subsection 3.1). By filtering the lightlike singularities with 
these properties in the above proposition we arrive to the following 

Proposition 8.3 A general spherically symmetric spacetime with metric |5U 
admitting, first, a natural number n > 1 such that lim( n _ > o.i?-s>o) d n m/du n (u, R) ^ 
(while, in case n / 1, limjy^-R-^o) d z m/du l (u, R) = for i = l,...,n — I) and, 
second, C n extensions rh and (3, will develop a persistent lightlike naked singularity 
as a consequence of gravitational collapse if e = +1, m(0, 0) = ; m,R (0, 0) < 1/2, 

[l-2m, fi (0,0)] 2 -16m, M (0, 0) > 0, 
d n rh , 

Ti— (0,0) > 0. 

Note that persistent lightlike naked singularities are also possible in the especial 
case fh,R (0, 0) = 1/2. The reader can write explicitly the conditions for them 



directly from figures 16, 17 and 18 



So far we have been working with a single local chart in U and assuming 
that certain extensions of m exist and possess a certain degree of differentiability. 
Therefore, the cases in which these requirements are not fulfilled have been excluded 
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up to now. However, in some cases, the method can be generalized to include 
different local charts or more general differentiability requirements. For example, 
particular models for collapsing stars with an initially regular center (rh(u < 0, R — 
0) = 0) in which the shell focussing eventually generates a persistent lightlike 
naked singularity (u = 0) could be constructed with the usual matching techniques 
(see [13j[26j for particular cases), which allow the avoidance of the differentiability 
requirements and the limitations of the use of a single local chart. It is only 
required that there is a regular spacetime modelling the initially regular interior, 
that there is a singular spacetime whose scalar invariant m satisfies the conditions 
in proposition |8.3 and that the two spacetimes can be matched through a non- 



spacelike hypersurface reaching R = at u = 0. 

In the next subsection we will also go beyond the limitations of dealing with 
a single local chart and with the required differentiability assumptions by working 
with a particular interesting novel application: Evaporating Black Holes that can 
develop persistent naked singularities. With these treatments our results could be 
used to cover most of the spherically symmetric solutions found in the literature. 
Nevertheless, let us remark that one can still find particular cases where even the 
generalized treatment would fail. An example of this would be the model in pi] 
which possesses a function m that is not even well defined at u — R — 0. 

8.3. Beyond the restrictions: Evaporating black holes developing persistent naked 
singularities 

Let us consider a future spacelike singularity that reaches a point, say at u = 0, 
where linv^^o m = and then it is followed by a regular R = -see ([3 Iff - for 



u > 0. If this happens, we will say that the singularity evaporates. In fact, this 
is the usually expected behaviour for the singularity of an evaporating black hole 



(EBH) [53] . (We have illustrated it with a particular complete model in figure 20). 

In this case there will be a future directed radial null geodesic starting at u = R = 

in the singularity. In this way, the singularity must be naked and this geodesic 

ff Note that the singularity disappears for u > provided that f3 also satisfies its own regularity 
conditions (pi) as it approaches R — for u > 0. 
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Figure 20. Typical Penrose's diagram for a collapsing star generating a BH that 
evaporates due to the emission of Hawking radiation. The grey zone corresponds 
with the interior of the star. The exterior of the star will usually contain radiation 
coming out from the star (the wavy arrows) and particles and radiation due to 
the Hawking effect. 



defines a Cauchy Horizon (CH) of the model. Sometimes this singularity is called 
an instantaneous naked singularity [29] since an observer crossing the CH could only 
detect a single flash of null radiation coming from the singularity. 

Now we would like to show that the requirement of regularity for u > in U 
allows the interesting possibility of the development of a persistent naked singularity 
at u — 0. Specifically the results in section [6] imply the following 

Proposition 8.4 A singular spherically symmetric spacetime with metric ([7]) 
admitting a natural number n > 2 such that lim( u _ > o,/?^o) d n m/du n (u, R) ^ 0, while 



lim (u ^ ,i?^o) d l m/du l (u, R) = for i 



1, ...,n 



1, admitting C n extensions m and 



(3 and with a R = future space-like singularity will develop a persistent lightlike 
naked singularity at R = u = if e = +1, 

m(0,0) =0, 

m,„(0,0) =0, 
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aVn (o,o)>o, 



du r 

provided that n is even, and one of these two options is satisfied 



a) 

dR K "'"' ~~ 2 



dm 1 . , 

(0,0) <- (25) 



dm 1 d % ( dm \ , . . 

r(0,0) = s ^ hnj 0,0)<0, (26) 



• b) 

dR^'" J 2 """ du l \dR ) 
where i is the lowest number such that the partial derivative is non-null provided 
that i < n/2. 

Note that the requirement n > 2 indicates that the generation of such light-like 
singularities is related to the presence of a scalar invariant m reaching its zero value 



slowly enough. In figure 21 we have sketched a procedure for constructing a future 
spacelike singularity developing a persistent lightlike naked singularity and followed 
by a regular center using the usual matching technique. On the other hand, in 



figure 22 we show a complete spacetime with a persistent lightlike globally naked 
singularity and a regular R = 0-center for u > 0. 
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Appendix A. Semi-Hyperbolic Critical Points 

Let us assume that the origin is an isolated critical point of a planar system. Then, 
if the linearization matrix A satisfies det A = and its trace T^O then the critical 
point will be semi-hyperbolic [TO]. The details on the theory regarding the qualitative 
behaviour of these points can be found in [9] [10] . We will summarize here only the 
main results: 
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a) 



o 
II 

Pi 

R = 
/wwwvw\k; 



b) 



Figure 21. The upper part of the graphic shows the second family of radial null 
geodesies for an m satisfying the requirements in proposition |8.4| together with 
its corresponding sketched conformal diagram. We have outlined two different 
lightlike hypersurfaces as our candidate for matching hypersurface: Ej and £p. 
In the first case, the only lightlike geodesic with a regular direction £j is chosen. 
The grey region is replaced with a region with a regular center and the result of 
the matching possessing an instantaneous naked singularity is shown in a). In 
the second case, one of the outgoing lightlike geodesies leaving u = R = with 
infinite slope is chosen as the matching hypersurface Sp. The darker grey region 
is replaced with a region with a regular center and the result of the matching 
possessing a persistent naked singularity is shown in b). 



The neighborhood of the origin can be divided into open regions called 
sectors which can be hyperbolic, parabolic or elliptic according to their topological 



equivalence with the figures Al-H, |A1| -P or Al-E, respectively. The trajectories 
which lie on the boundary of a hyperbolic sector are called separatrices. 

Given a system possessing an isolated semi-hyperbolic critical point there is a 
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Figure 22. Penrose's diagram for a collapsing star generating a BH that 
evaporates due to the emission of Hawking radiation. The scalar invariant m 



satisfies the requirements in proposition 8.4 for u < and a regular center is 
proposed for u > 0. The resulting diagram contains a persistent globally naked 
singularity: The light-like R = u = singularity. 




Figure Al. (H) Hyperbolic sector, (P) parabolic sector and (E) elliptic sector. 



suitable linear transformation that allows us to write it as: 

f = P2{x,y) 

%=y + Q 2 ( X} y) } 

where P 2 and Q 2 are functions satisfying P 2 (0, 0) = Q 2(0,0) = P 2 , x (0, 0) = 
Q2,x(0, 0) = P2, y (0,0) = Q2, y (0,0) = 0. By the implicit function theorem, the 
equation y + Q2(x,y) = has a solution y = <p(x) in a neighborhood of O. If we 
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define the function if)(x) = -P 2 (x, (p(x)) then its truncated series expansion will have 
the form 



where n > 2 and A n ^ 0. Then it can be shown ([9], Theor.65, p. 340) that 

• If n is odd and A n > 0, (0, 0) is a topological node, i.e., there is a trivial sectorial 
decomposition consisting in only one parabolic sector. 

• If n is odd and A n < 0, (0, 0) is a topological saddle, i.e., four hyperbolic sectors 
separated by four separatrices. Two of these separatrices tend to (0, 0) in the 
directions and ir, the other two in the directions 7r/2 and 37r/2. 

• If n is even then (0,0) is a saddle node, i.e., one parabolic and two hyperbolic 
sectors separated by three separatrices. 

— If A n < 0, the hyperbolic sectors contain a segment of the positive x-axis. 

— If A n > 0, the hyperbolic sectors contain a segment of the negative x-axis. 

Appendix B. Nilpotent critical points 

Let us assume that the origin is an isolated critical point of a planar system and 
that A (its linearization matrix) is not the zero matrix, but det A = and its trace 
satisfies T = then the critical point will be nilpotent [10J. Again we refer the reader 
to [9] pH] for details on the theory of the qualitative behaviour of these points. We 
will summarize here only the main results: 

Given a system possessing a nilpotent critical point there is a suitable linear 
transformation that allows us to write it as: 



where P 2 and Q2 are functions satisfying P 2 (0,0) = (52(0,0) = P2, x (0, 0) = 
<52>z(0,0) = P^jy (0,0) = Q^jy (0,0) = 0. By the implicit function theorem, the 
equation y + P2(x,y) = has a solution y = (p(x) in a neighborhood of O. If we 
define the function ip(x) = <5 2 (a;, <p(x)) then its truncated series expansion will have 
the form 



i>(x) = A n x n + ... 



(A.l) 




ip(x) = a k x k + ... 
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where k > 2 and 7^ 0. On the other hand, if we define a(x) = P2, x (x,<p(x)) + 
Q 2 , y (x, <p(x)) then either its truncated series expansion can be written as 



where b n ^ or a(x) = in which case b n = Vn. Then it can be shown that 

• If k is odd (j9], Theor.66, p. 357) then we define p = (k — l)/2 and A = 
# + 2(fc + l)a fc 



— If dk > then O is a topological saddle point. 

— If <2fc < then O is a 

* Focus or center if either 1) b n = 0, 2) 6 n ^ and n > p or 3) 6 n 7^ 0, 
n = p and A < 0. 

* Topological node if either 1) n is even, 6 n 7^ and n < p or 2) n is 
even, b n ^ 0, n = p and A > 0. 

* Equilibrium state with an elliptic region, i.e., an elliptic sector and a 
hyperbolic sector separated by two separatrices, if either 1) n is odd, 
b n 7^ and n < p or 2) n is odd, 6 n 7^ and n = p and A > 0. 

• If & is even ([9], Theor.67, p. 362) then we define p = k/2 and the critical point 
is a 

— Cusp (i.e., two hyperbolic sectors separated by two separatrices) if either 
b n = or b n 7^ and n > p 

— Saddle node if 6 n 7^ and n < p 
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